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Abstract 

Fission chambers were first made fifty years ago for neutron detection. At the 
moment, the French Atomic Energy Commission (CEA-Cadarache) is developing a 
sub-miniature fission chamber technology with a diameter of 1.5 mm working in the 
current mode [1]. To be able to measure intense fluxes, it is necessary to adjust the 
chamber geometry and the gas pressure before testing it under real neutron flux. In 
the present paper, we describe a theoretical method to foresee the cur rent- volt age 
characteristics (sensitivity and saturation plateau) of a fission chamber whose ge- 
ometrical features are given, taking into account the neutron flux to be measured 
(spectrum and intensity). The proposed theoretical model describes electric field 
distortion resulting from charge collection effect. A computer code has been devel- 
oped on this model basis. Its application to 3 kinds of fission chambers indicates 
excellent agreement between theoretical model and measured characteristics. 

Key words: fission chamber, electric field distortion, fission fragment, ionization, 
ion-pair 



1 INTRODUCTION 



A fission chamber is a device that measures a neutron flux and that is weakly 
disturbed by the gamma flux [2, p 503]. 

In fig. 1, we show the basic principle behind the use of a fission chamber as 
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Fig. 1. A schematic view of a fission chamber. 



a neutron detector. A thin layer of fissile material ( 235 U for instance) is de- 
posited on the anode. When a neutron traverses the detector, it may collide 
with a fissile atom (the probability is proportional to the fission reaction cross 
section) and the resulting reaction generates two charged particles (Fission 
Fragment : FF) produced by fission. Those charged particles cross the argon 
atmosphere between the two electrodes of the chamber (anode and cathode) 
and ionize the gas on their trajectory creating an average of 0.3 ion-pairs/nm 
(Ar + , e~). A DC voltage of a few hundred volts being applied between the 
two electrodes, those electrons and ions drift across the gas in response to the 
applied electric field producing a signal that can be amplified and processed. 
The voltage-current characteristic of a fission chamber can present a satura- 
tion plateau. It means that, in a voltage interval, the current delivered by 
the chamber is roughly constant. The saturation current is proportional to 
the neutron flux, that is the reason why we need to know when the cham- 
ber works in saturation conditions. But, for a given geometry, the so called 
voltage interval depends on the neutron flux. The main aim of this study is 
to determine the relevant parameters of the chamber and the optimal voltage 
that answer the neutron flux requirement of the user. 

Firstly, we will study quantitatively the conversion of neutrons into fission frag- 
ments and into ion-pairs. Secondly, we will spend some time to understand 
the charge collection and the creation of the signal. Then, we will compare the 
results predicted by our model with those obtained experimentally. 



2 ION-PAIR CREATION 



2 



2.1 Neutrons make Fissions 



We are dealing with the fission reaction that turns an uncharged particle 
(neutron) into charged particles (FF) which are easier to detect : 



Let fi s be the surface mass of fissile product (for a subminiature fission chamber 
fi s = 760 fig ■ cm -2 ). This parameter decreases during the life of the chamber, 
because fissile atoms are not regenerated. This is our stock of fission fragments. 
The rate of fission reactions per surface unit is given by 



(the subscript fst stands for fission per unit surface per unit time). My is 
the molar mass of 235 U, N A is the number of Avogadro, er(e) is the fission 
cross section for neutron of energy e and ^(e) is the energetic density of the 
incoming neutron flux. For the next steps we will call x — Io° (T ( e ) c ' ) e( e ) de. A 
same value of x can be generated by various neutron spectra. We just have to 
find the right spectra of incoming neutron ^(e). This could cause a problem 
as we are going to see in the next part. 

2.2 Fissions make fission fragments 

For example : take two different spectra of incoming neutrons that generate 
the same fission rate x- Let us consider one flux of incoming neutrons of energy 
0.04 eV and one of energy 1 MeV. For 235 U, 



Then a flux of 1.0 10 12 n • cm -2 • s -1 for 0.04 eV neutrons and a flux of about 
5.7 10 15 n ■ cm -2 ■ s -1 for 1 MeV neutrons result in the same reaction rate x '■ 
X = 5.7 10 _10 s _1 and the same number of fission fragments is created. The 
problem is that, a priori, the nature and the energy of these fragments is 
not the same as the ionization rate in the gas. Fortunately, the distribution 
of fission fragments (fig. 2) do not vary significantly with the energy of the 
incoming neutrons. Moreover, since the average energy of the 2 FF is about 
160 MeV, we can say that even a difference of 15 MeV in the energy of incoming 
neutron induces an error of less than 10 %. So, we will suppose that the fission 
fragments distribution and the fission fragments energy distribution do not 
vary significantly for incoming neutron energies between MeV and 15 MeV. 




(2FF, 2.3neutrons, energy) 




(j(0.04 eV) = 570 barn and tr(l MeV) = 0.1 barn. 
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Fig. 2. distribution of 2ci5 U fission fragments. 
2. 3 Fission fragments make ion-pairs 

One fission gives birth to two fission fragments. Their average energy is around 
90 MeV for the lightest and 60 MeV for the heaviest. If the energy of the 
incident neutron is negligible compared to 90 MeV and 60 MeV it is accurate 
to say those FF emerge in opposite directions. Thus only one FF participates 
in the gas ionization, the other one being absorbed by the anode. 



2.3.1 Calculation of the number of ion-pairs created per unit time per unit 
volume 

We want to calculate the density of ion-pairs created per unit time (N(r)) at a 
distance r from the anode axis (see on fig. 3). We assume the hypothesis that 
FFs are emitted in every direction with the same probability. Let N ft be the 
number of fission reactions on the surface element dS per unit time. 

N ft = N fst dzr a dfl 

Let l c (x) be the average number of ion-pairs created by a fission fragment per 
unit of length travelled per unit pressure, x representes the distance travelled 
from the creation point. We can calculate the number of ion-pairs created per 
unit of time in the volume element dr = r'd^dr'dz around the point P. This 
number is 

dN = N ft cos(^)^pl c (r')dr' 
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where p is the pressure of the chamber. Thus, the density of ion-pairs created 
in P because of dS is 

dN N fst r a 

— = -cos(V>)d0 pl c (r ) 

ar n r 

In the case of a sub-miniature fission chamber the trajectory of fission frag- 
ments is short enough to consider l c (r) to be constant equal to Iq Thus we can 
say that the density of ion-pairs created at distance r from the axis per unit 
time is 




Fig. 3. Schematic cross section of the Anode. 



2.3.2 Justification of\ c {r) constant equal to \q 

The results of this part are due to the program SRIM of J.F Ziegler which 
calculates the stopping and range of ions in solids, liquids and gases [3]. In 
fig. 4, we show the energy loss of two typical fission fragments : 95 Mo at 
90 MeV and 139 La at 60 MeV. The space separating the electrodes of all types 
of fission chambers is no longer than a few milimeters. For such a distance, 
we can notice that the energy loss by ionization in argon is roughly constant 
equal to 4.5 keV.yum -1 and so is l c (r). Suppose that the creation of one ion- 
pair requires 23 eV (ionization energy of argon). We can say that the number 
of ion-pairs created by one typical FF per unit of length through the gas is 
Iq = 2. 10 5 ion-pairs/mm/bar. This hypothesis will be the first source of error. 
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Fig. 4. Energy loss for Mo and La in 1.0 bar Argon. 
2.3.3 Asymptotic hypothesis 

It would be quite difficult to work with (1). It simplifies singularity if, in- 
stead, we take the asymptotic version (r — > r a ) or (r — ► +oo). Suppose 
(r c — r a ) <C r a then r remains in the vicinity of r a and we can assume that 
N(r) is constant and equal to its value in r a 

N(r) = N = ^N fst pl 

Suppose, now, (r c — r a ) r a then we can take with a good accuracy the equiv- 
alent of N(r) at infinite, 

N(r) = N fst pl^ = t 
r r 

In fact, the first approximation fits with the case of plane electrodes, and the 
second squares with cylindrical geometry of the chamber. Those hypothesis 
will be the second source of error. 



3 PHYSICAL PHENOMENA IN THE GAS 



Let N(r) be the number of ion-pairs created by fission fragments per unit of 
time and volume in the gas (it depends on x — I °"(e) ( t | e( e )de as seen before). 
We will suppose this number to be equal to N in the entire volume of the gas, 
because most fission chambers verify the geometrical hypothesis (r c — r a ) <C r a . 
This approximation will give us only an order of magnitude. But, from an 
experimental point of view, it would be utopic to expect better than this. 
Indeed, the geometry of a fission chamber varies from one to the other. It is 
impossible to get a perfectly cylindrical geometry, and this may cause some 
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difference of charge collection between two similar models of fission chambers. 
Various physical phenomena can occur in the gas [4]. Comparing experiment 
and theory, we noticed that only certain effects could be taken into account. 



3.1 Primary ionization 

The primary ionization corresponds to the ionization reaction 
(FF, Ar) — > (FF, Ar + , e - ). The fission fragments ionize the argon and generate 
ion pairs (Ar + ,e~) and they lose energy by this process (see on fig. 4). 



3.2 Electrical drift and temperature 

The DC voltage applied between the two electrodes generates an electric field 
that is responsible for the movement of electrons towards the anode and ions 
towards the cathode. This is the electrical drift. The respective velocity of 
electrons and ions in such a field are //_ E and ji + E where /z_ and ji + are the 
mobility of those particles in the gas . In a gas, the mobility of a charged 
particle is assumed to be inversely proportional to the pressure of the gas [4]. 
For argon : 



Since the pressure is proportional to the temperature, we can say that the 
mobility varies as the inverse of the temperature. This will be the only effect 
of the temperature modelized. 



3.3 Recombination 

Since Ar + and e~ are oppositely charged particles, it is also possible that a cer- 
tain amount of those particles recombine into mother particles (e~, Ar + ) — > Ar. 
This is the recombination reaction [5] . The rate of recombination depends on 
the local density of ions and electrons and is given by 



H_ = O.lSm^bar.V^.s- 1 



= o.mcrW.barAr 1 ^ 1 




The tabulated value of a is assumed to be 10 12 SI for argon 
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3.4 Threshold condition 



Now, suppose that the electric field is strong enough to give the free electrons 
a kinetic energy of a few eV (order of magnitude). The secondary ionization 
reaction (e~,Ar) — > (2e~,Ar + ) is then possible. The thermal kinetic energy 
is about 0.04 eV. It means that the kinetic energy of a few eV of electrons is 
mainly due to the electric field : Kinetic energy = |m e (/i_Es/p) 2 is about a few 
eV. For a 1.0 bar pressure of Argon we can calculate the electric field threshold 
(E$) that permits the secondary ionization reaction : 

Es = some units of 10 6 V • rrT 1 

In this way, this hypothesis can explain the fact that the end voltage of the 
saturation plateau is higher when the pressure is higher (ref. part 2.1). Indeed, 
in order to maintain the kinetic energy constant, we must keep the ratio E$/p 
constant. It means that Es is proportional to the pressure of argon. In our 
model, we took, for a 1.0 bar chamber, E s = 1.7 10 6 V • m -1 . 



4 CHARGE COLLECTION 

4-1 Minimum voltage for steady state solution 

The relevant signal of a fission chamber in the current mode is the RMS 
value of the current. We are not interested in the noise generated by the 
chamber. Thus, we are going to study the stationary state of the chamber. If 
we want a steady state solution to exist, all phenomena must balance. Suppose 
that the only effects we take into account are the primary ionization ( this 
is the dominant effect on the saturation plateau ) and the electrical drift. 
As soon as recombination is of the order of magnitude of primary ionization 
this hypothesis will not be valid anymore. These hypotheses will enable us 
to compute easily the electron density, the ion density and the electric field. 
We expect a minimum value of the voltage under which the electric field is 
not strong enough to evacuate the charges created by the primary ionization. 
Thus, below this minimum voltage a steady state solution could not exist. We 
will evaluate the validity domain of these hypotheses by comparing the value 
of «n_n + a posteriori and the value of N. But, for now, let us compute the 
steady state solution of this problem with these hypotheses 

-^d r (r/i_En_) = N(r) = N (2) 
^d r (r /U+ En + ) = N(r) = N (3) 
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1 e 
-d r (rE) = — (n + - n_) 

f Co 



(4) 



Since the voltage is maintained constant we have an integral condition on the 
electric field whose structure depends on the space charge distribution (4) 



E(r)dr = V p (5) 



We need two boundary conditions on n + and n_. Ar + drift towards the cathode 
and e~ towards the anode. For a steady state solution, in order to avoid an 
accumulation of Ar + at the cathode and e~ at the anode, n_ and n + must verify 

n + (r a ) = (6) 
n-(r c ) = (7) 



The system of non linear differential equations (2) to (4) with conditions (5) 
to (7) is self-sufficient. They induce 



Thus, thanks to (4) 

I dr((rEf) = No i((^il) + (^)) 

r e V fi + /I. J 

We can integrate this differential equation and then we get 

E(r , No , C ) H N.f(^ + ^) + ^ 
4e V fi+ V- ' H 

where C is a constant of integration whose value can be calculated with the 
condition (5). We notice that we cannot find a C for every V p . This is what 
we expected. For a given geometry of the chamber and for a given value of 
Nq, there is a minimum voltage V m j n (N ) under which there is not any steady 
state solution for the problem. The value of C = C m ] n for this minimum voltage 
is given by the condition E(r a ) = 0. For every C > C m ; n , the electric field is 
nonzero, which is a necessary condition for n_ and n + to exist (see (9) and 
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(10)). 




) 



Thus we can say that the begining voltage of the saturation plateau is greater 
or equal than 



whatever a is (the recombination constant) we cannot find a saturation below 
this limit. We can notice the importance of space charge effect at the begining 
of the saturation plateau. 

4-2 The saturation domain on a voltage-flux graph 
4-2.1 Recombination limit : V rec 

The saturation current corresponds to the collection of all electrons created 
by the primary ionisation. So, l sat is equal to e N(r)27r r L dr, that is to say 



We can notice that l sat depends on the product (r c — r a )(r c + r a ) whereas for 
a chamber for which (r c — r a ) 3> r a , N(r) would be (5/r and thus l sat would 
depend only on (r c — r a ). The limitations of the saturation plateau are due 
to the recombination and the threshold condition. The recombination modifies 
the saturation current by the following quantity 



This quantity can be evaluated by replacing n + and n_ by their value calculated 
in (4.1). Since the value of a is very weak we can consider that S\ rec becomes 
important, relatively to l sat , at the vicinity of V m j n (N ). Indeed, at this point 
E(r a ) = 0, so n + (r a ) and n_(r a ) tend to infinity (theoretically speaking). So 
is 5\ rec . That is why we assume that the beginning voltage of the saturation 
plateau V rec is equal to V min (N ). 

4-2.2 Threshold limit : V max 

Suppose that over a certain electric field threshold (E$) there is an electronic 
multiplication that ends the saturation plateau. We want to determine, for a 
given value of No, the minimum voltage corresponding to the beginning of this 
phenomenon. It appears when either E(r a ) or E(r c ) becomes greater than E$ 




l sat = 7reLN (r c 2 -r a 2 ) 



(10) 




v 
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since the function E(r) is concave. 

Suppose N = 0. There is no ion-pair created in the volume of the gas. So there 
is no space charge problem. Thus, the electric field has the well known form 
in 1/r (see graph. 1 on fig. 5). We reach the threshold condition if E(r a ) = E s 
that is to say Q = E|r a and then V max (0) = / E(r)dr = Esr a / dr/r. 
Suppose now that we increase the value of No- Since ion-pairs are created in 
the gas, space charge may cause electric field distortion. On fig. 5 , we can see 
the evolution of the electric field with N (the voltage is maintained constant). 
We can meet three cases at the threshold condition : either E(r a ) > E(r c ) and 
then E(r a ) = E s (graph. 1 on fig. 5), E(r a ) = E(r c ) and then E(r a ) = E(r c ) = E s 
(graph. 2 on fig. 5), or E(r a ) < E(r c ) and then E(r c ) = E$ (graph. 3 on fig. 5). 
The first case involves that 

2.2 , e ,/r* , (2r 2 -r 2 ) 



d = Egr +— N r ^ + 

4e \fj,+ li- 

and the corresponding voltage 

V maxl (N )= J E(r,N ,d)dr 
The second and third case involve that 

r r2 2 e m ~ 2r a) r c 

C 2 = E 2 r c 2 - —N r c 2 ^ ^ - 

4e V /!+ li- 

and the corresponding voltage 

V max2 (N )= J E(r,N ,C 2 )dr 

We can notice, because of the expression of Q and C 2 , that space charge 
becomes important for high value of Nq. Moreover this space charge effect 
is responsible for the maximum value of N 0max over which saturation cannot 
exist. 



4-3 The use of these results 



Suppose we want to know the saturation domain of a given fission chamber 
for a given incoming neutron spectra. We have to calculate x — I (7 ( e ) c ' ) e( e )de. 
For this purpose we can use a Monte Carlo code, like MCNP, developed by 
the Los Alamos National Laboratory, to simulate the collision of neutrons 
with the fission chamber and the fissile product inside, and then to calculate 
X and N . Then, knowing the geometry of the chamber, we can get V rec (N ), 
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Fig. 5. Electric field for various values of No at constant voltage. 

which is the begining voltage of the saturation plateau, and V max (No), which 
is the end of the saturation plateau. We can sum these results up in a graph. 
The abscissa is the voltage and the ordinate is either Nq or the current or 
the thermal flux of neutrons, depends on what we want to know. On fig. 6 
and fig. 7 we draw various saturation domains, each being limitated by two 
curves. These curves represent V rec and V max . We wrote a very simple code 
that permits us to generate the graph that represents V min (N ) and V max (N ). 
The inputs of this code are : r a , r c , L (geometry) and M f (fissile product). 
The outputs are three arrays of points : V m i n (N ), V max (No) and N . 



5 MODEL VERSUS TECHNICAL DATAS 



We want to compare the results we get by the use of our model with the 
characteristic of two industrial fission chambers manufactured by PHOTON IS, 
imaging sensors : CFUR 43 and CFUF 43. Let us begin by the first one. The 
table 1 gives us the geometrical and physical features of the chamber. 
We put r c , L, pressure and \i a to their actual value in our code. PHOTONIS, 
imaging sensors writes in the CFUR43 manual that the optimum polarisation 
is 150 V. Let us say that we do not know r a . We want to find its value with 
the help of our code. We draw on the following graph (fig. 6) various domain 
of saturation of CFUR 43 for various value of r c — r a (inter-electrode gap). 
We can remark that values of r c — r a between 0.2 mm and 0.3 mm fit with 
the requirement of a 150 V optimum polarisation. That is exactly^ the range 
value given by PHOTONIS, imaging sensors 



2 the exact value of r a is an industrial secret of PHOTONIS 
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Fig. 6. Variation of the saturation domain with the inter-electrode gap for CFUR 
chamber 

You could say that it only works for CFUR 43. So let us see it works with 
another chamber with different size : CFUF 43 (see table 2). The result given 
by our model can be visualised on fig. 7. 



Saturation domain for CFUF 43 for various 
inter-electrode gaps 




Supplied voltage (V) 



Fig. 7. Variation of the saturation domain with the inter-electrode gap for CFUF 
chamber. 

We can also compare the sensitivity of these fission chambers (given by PHO- 
TONIS) with the one calculated with our model. The sensitivity is the ratio 
of the current delivered by the chamber and the thermal neutron flux 

The various results are summarized in table 3. Tabulated and calculated values 
are in very good agreement. We can deduce from these two figures that an 
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increase of (r c — r a ) induces a decrease of the maximum flux detectable. It also 
induces an increase of the value of the optimum voltage for saturation, and of 
the width of the saturation plateau. 



6 MODEL VERSUS SILOE RESULTS 



In this part we present the 1.5 mm diameter subminiature fission chamber 
manufactured at the CEA-Cadarache, France (see table 4). We tested two home 
made 1.5 mm diameter chambers : a 1.1 bar and a 4.0 bar chamber (this is the 
pressure of argon in the chamber). These experiments were carried out at 
the CEA/Grenoble (France) at the SILOE facility (experimental nuclear power 
plant). The aim of these experiments was, first, to make some general empirical 
observations, and second, to collect enough data to validate our model. 



6.1 Comparison between a 1.1 bar and a 4.0 bar fission chamber 



We compared the characteristics of both fission chambers in the same thermal 
flux of about 3.4 10 12 n.cm _2 .s _1 . We can see in fig. 8, as expected, that the 
saturation current for a 4.0 bar fission chamber is about four times greater 
than the one for a 1.1 bar chamber. Actually, the accurate ratio is 3.5. But 
we can also notice that the saturation plateau, for a 1.1 bar chamber, begins 
at a voltage lower than the one for a 4.0 bar chamber. Thus, we can deduce 



Neutron flux * lrermol =3.4 10 ' n s 1 cm 



3 

>n 80 
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Fig. 8. Characteristic of two chambers for the same neutron flux. 

from those remarks the following rule : if, for a given neutron flux, the 
voltage to reach the saturation plateau is too high, we can reduce 
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the argon pressure in the chamber, but as a consequence, we also 
reduce the saturation current. This rule is verified with our model. On fig. 
9 and fig. 10 one can see the superposition of experimental points and of V rec 
(line representing the beginning of the saturation domain). On these figures, 
"?" means that we didn't know the thermal neutron flux at this location in the 
core. Once more, the model is in good agreement with the experimental data. 
It predicts with a good accuracy the beginning of the saturation domain. 



'.1 bar fission chamber 
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Fig. 9. Characteristic of a 1.0 bar chamber for various neutron fluxes. 



4.0 bar fission chamber 
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Fig. 10. Characteristic of a 4.0 bar chamber for various neutron fluxes. 

For a 1.1 bar chamber we measured characteristics up to 270 V. Above this 
voltage the current intensity increased very quickly. That is why we did not 
measure further in order not to destroy the chamber. This is the same reason 
why we did not carry out the experiment above 430 V for the 4.0 bar cham- 
ber. We calculated the sensitivity of each chamber with measured data (for 
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4> th = 3.4 10 12 n.cm .s , the l.lbar chamber gives l sa t — 37/xA and the 4.0bar 
chamber gives l sat = 120/iA) and with our model. The results are summarized 
in table 5. The order of magnitude is the same. The discrepancy between mea- 
sured and calculated values is consistent with the uncertainty of the pressure 
and geometry. 



7 CONCLUSION 



The model described in this paper is in good agreement with technical data for 
existing fission chambers CFUR 43 and CFUF 43, and also with experimental 
data obtained with home made subminiature fission chambers. The model also 
takes into account the temperature, which is necessary for the use of a fission 
chamber in the core of a nuclear power plant. Finally, this study revealed 
two important behaviors of fission chamber : firstly, in order to reduce the 
value of the minimum voltage of saturation, we can reduce the pressure in 
the chamber, but, as a consequence, we reduce the sensitivity of the chamber. 
Secondly, we learned that increasing the inter-electrode gap leads to a decrease 
of the maximum flux detectable and an increase of the optimum voltage of 
saturation and of the width of the saturation plateau. 
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Table 1 

Geometrical and physical features of CFUR 43. 



CFUR 43 Features 



r c 1.25 mm 

L 13 mm 



pressure 1.1 bar 

Us 200 ng ■ cm~ 2 



Table 2 

Geometrical and physical features of CFUF 43. 



CFUF 43 Features 



r c 1.78 mm 

L 27 mm 



pressure 1.1 bar 

fi s 160 fig ■ cm~ 2 



Table 3 

Comparaison of tabulated and calculated sensitivity for CFUR 43 and CFUF 43. 



Sensitivity of fission chambers 




CFUR43 


CFUF43 


Stabulated(A/n.Cm 2 .s x ) 


3. 10- 18 


1. lO" 17 


Scalculated(A/n.Cm~ 2 .S- 1 ) 


3.3 10~ 18 


0.9 10~ 17 



Table 4 

Geometrical and physical features of subminiature fission chamber. 



Subminiature Fission Chamber Features 


anode radius r a 


0.35 mm 


cathode radius r c 


0.65 mm 


length of sensitive part L 


12 mm 


pressure of argon 


1.1 bar + 4.0 bar 


surfacic mass of fissile product ^x s 


760 jUg • cm~ 2 
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Table 5 

Comparaison of measured and calculated sensitivities for 1.1 bar and 4.0 bar sub- 
miniature fission chambers. 



Sensitivity of subminiature fission chambers 




1.1 bar chamber 


4.0 bar chamber 


S me asured(A/n.Cm 2 .S X ) 


1.1 1(T 17 


3.6 10~ 17 


S C alculated(A/n.Cm" 2 .S- 1 ) 


0.7 1(T 17 


2.8 10~ 17 
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